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ABSTRACT 

Globular clusters were thought to be simple stellar populations, but recent photo- 
metric and spectroscopic evidence suggests that the clusters' early formation history 
was more complicated. In particular, clusters show star-to-star abundance variations, 
and multiple sequences in their colour-magnitude diagrams. These effects seem to be 
restricted to globular clusters, and are not found in open clusters or the field. In 
this paper, we combine the two competing models for these multiple populations and 
include a consideration of the effects of stellar collisions. Collisions are one of the 
few phenomena which occur solely in dense stellar environments like (proto-)globular 
clusters. We find that runaway collisions between massive stars can produce mate- 
rial which has abundances comparable to the observed second generations, but that 
very little total mass is produced by this channel. We then add the contributions of 
rapidly-rotating massive stars (under the assumption that massive stars are spun up 
by collisions and interactions), and the contribution of asymptotic giant branch stars. 
We find that collisions can help produce the extreme abundances which are seen in 
some clusters. However, the total amount of material produced in these generations 
is still too small (by at least a factor of 10) to match the observations. We conclude 
with a discussion of the additional effects which probably need to be considered to 
solve this particular problem. 

Key words: globular clusters: general - stellar dynamics 



1 INTRODUCTION 

Globular clusters have long been viewed as the epitome of 
simple stellar populations. Their stars have a common age, 
a common distance, and a common metallicity; there is no 
interstellar gas and little else to get in the way of studying 
the stars directly. These systems are the closest we can come 
to a "controlled experiment" in stellar astrophysics, and as 
such they have proved incredibly valuable for studies of both 
stellar evolution and stellar dynamics in the past. 

However, in recent years, cracks have been appearing 
in this simple picture. Both photometric and spectroscopic 
studies of clusters have started to unearth puzzles and in- 
consistencies. At first, these problems were thought to be an 
oddity in one particular cluster, or evidence of some details 
of stellar physics that we didn't quite understand. Over the 
last five years or so, however, it is becoming clear that our 
picture of a globular cluster needs to change. They cannot 
have formed instantly out of a single molecular cloud, re- 
moving all their leftover gas immediately, and then evolved 
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passively for the next 10 billion years. Their history is more 
complicated. 



1.1 Observational Background 

Hints that something strange was going on in globular clus- 
ters came first from spectrosco pic studies of their red giants. 
For a n excellent review, see IGratton. Sneden fc Carrettal 
The general results have not changed since that re- 
view was written, although we now have observations of 
more stars per cluster, and more stars observed with high- 
resolution spectra. Most globular clusters have a constant 
iron and iron-peak ele ment abundances, with the notable ex- 
ception of w Centauri (|Freeman fc Rodg ers 1975 ) and hints 
of a very small spr ead in M22 and M92 (jMarino et alj|2009l ; 
lLanger et al.|[l99a ). However, it has been known since the 
late 1970s that lighter metals, particularly carbon, nitro- 
gen and oxyg en, do vary from star to star in many clusters 
Other light elements also show star-to-star 
variations in clusters, including Na, Al and Mg. The most 
striking piece of evidence to date that this abundance vari- 
ation phenomenon occurs in all clusters is the spectroscopic 
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study of sodium and oxygen bv lCarretta et alj (|2009al ibh. A 
general anti-correlation is seen, with oxygen-depleted stars 
having higher sodium abundances. Aluminum and magne- 
sium have been studied in fewer clusters, but they als o show 
similar trends (IShetroneHl99rj ; ICarretta et al, 2009b), with 
a large range in aluminum, a smaller spread in magnesium, 
and some evidence for high-aluminum stars having lower 
magnesium abundances. 

ICarretta et al.l |2008l ) also looked for correlations be- 
tween the extent of the Na-0 anti-correlation and cluster 
properties. The strongest correlation was between the extent 
of the correlation and the maximum temperature of stars on 
the horizontal branch, in the sense that clusters with very 
hot (or blue) horizontal branches had a large spread in Na 
and O. They also found a weaker trend with the total mass 
(or magnitude) of the cluster, and a trend with galactic or- 
bit, in the sense that clusters which have spent more of their 
lifetime in the outskirts of the halo and not interacting with 
the galaxy, have a larger spread in abundances. 

One of the main reasons that globular clusters have 
been so useful to stellar astrophysics is that their colour- 
magnitude diagrams are very clean. Other than the blue 
straggler stars and a few other un usual objects, the star s 
fall onto a single isochrone (see e.g. ISaraiedini et al.ll2007l ). 
Again, the exception has been uj Centauri. Its colour- 
magnitude diagram shows a large amount of structure be- 
yond a single age/composition isochrone. The spread of the 
giant branch of this cluster was understood in the context 
of the measured differences in iron abundance of the stars, 
and so it did not come as much of a surprise to the com- 
munity when the main sequence of this c luster was found t o 
contain more than one sequence as well (|Bedin et al.| [2004). 
However, the most intriguing result from this detailed study 
of lj Centauri was the determination that the iron abun- 
dance of the bluest main sequence was not, as one would 
expect, the lowest in the cluster, but was in fact the highest 
l|Piotto et al.ll2005l ). The only possible way to reconcile the 
spectroscopic abundances with the photometric information 
was to infer a high helium content for these stars, perhaps 
even as high as Y=0.4. Because the abundance anomalies 
are found in light elements only but not iron, researchers 
have been casting this as a problem of 'pollution' from an 
early generation of stars in the cluster. 

The globular cluster community finally came to re- 
alize that mul t iple populations were ubiquitous when 
iD'Antona et all (|2005h discovered an intrinsic spread in the 
main sequence NGC 2808, which was later confirm ed to 
be three separate main sequences (|Piotto et al.l 12007). The 
turnoff region in this cluster is quite tight without much evi- 
dence for a spread, but starting about one magnitude below 
the turnoff, the sequences become obviously separated. 

At about the same time, we saw the first CMDs show- 
ing multiple subgiant branches (e.g. the ACS observation of 
NGC 1851 (|Milone et alj 120081 )'). Similar obse rvations have 
since been seen in M22, NGC 6388, and M54 (|Piottoll2009l ) 
and many intermediate age clusters in t he LMC and SMC 
l|Mackev et a.1.1 120081 : iMilone et alj|2009l ). NGC 1851 is an 
interesting case, is that it does not show any evidence for a 
splitting of the main sequence, even in very careful observa- 
tions including proper-motion cleaning of the CMD. It does, 
however, show a split in its giant branch when observed in 



U - I (|Han et alj|2009l) which is not seen in V - I, high- 
lighting the importance of observing in multiple bands. 

Understanding the horizontal branch (HB) morphology 
in globular clusters has been a problem for decades, com- 
monly referred to as the 'second parameter problem.' Be- 
cause the position of HB stars in the CMD is sensitive to 
both helium and metal abundances, they are a very useful 
population with whi ch to discuss both pollut ion and multi- 
ple populations (e.g. ID'Antona fc Caloill2008l ). 

1.2 Possible Explanations 

Many of the models to date have attempted primar- 
ily to_jmdfirst<m(i^^ The work 
of iPrantzos. Charbonnel."fc Iliadisl (|2007l) showed that the 
general abundance patterns could be explained by invok- 
ing the products of hot hydrogen burning, at a temperature 
of 70-80 x 10 6 K. At these temperatures, the hydrogen is 
burned in a series of cycles - the standard CNO cycle, as 
well as the neon-sodium cycle and the aluminum-magnesium 
cycle. Conveniently, the second two cycles work to give the 
abundance patterns that are needed: higher sodium, lower 
Mg and higher Al. This processed material must be removed 
from the star before helium burning can proceed. Otherwise, 
the helium will be converted to carbon and oxygen. We will 
not have the large amount of helium that is required at the 
surface, and the constraint of (almost) constant C+N+O 
abundance will be violated. 

A number of polluters have been proposed. The lead- 
ing contender is a population of intermediate-mass (3-10 
Mp ) asymptotic giant branch (AGB) stars (for a review, 
see iRenzinil 120081 ). These stars can reach the appropriate 
high temperatures in their hydrogen burning shells, and the 
processed material is brought to the surface by the outer 
convection zone as it reaches into the burning shell during 
the hydrogen burning portion of the thermal pulse phase. 
AGB stars have strong but low- velocity winds, which means 
that any mass that is removed from the star has a good 
chance of remaining in the cluster. Indeed, there is recent 
evidence that giant stars in this mass range do show surface 
abun dances which are consistent with th e necessary pollu- 
tion (|Villanova. Carraro. fc Savian 

An alternative source of pollution is a pop ulation of 
fast-rotating massive stars l|Decressin et alj 2007). Massive 
stars reach sufficiently high temperatures in their hydrogen- 
burning cores. These stars also have substantial winds, but 
under normal circumstances, the core regions are not ex- 
posed until very late in the stars' lives, well past the helium- 
burning phase. If a star is rotating rapidly, however, then 
meridional circulation and other rotational instabilities will 
mix material from the core to the surface, bringing these 
hydrogen-burning products up to a region where the wind 
can take them away from the star. Also, if the star is rotat- 
ing rapidly enough that it is close to its break-up velocity, 
material can escape from the equator of the star in a slow 
outflowing disc or slow wind. This material in particular has 
a slow enough velocity that it will stay in the cluster, and 
there are even indications that low mass stars could form in 
situ in the disc. 

Both AGB and fast-rotating massive star models have 
some difficulties in explaining the pollution of globular clus- 
ter gas and the formation of the second generation of stars. 
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In both cases, the amount of material that is required to 
form the second generation is quite large. Observations of 
clusters such as NGC 2808 suggest that the mass of the sec- 
ond generation is approximately equal to the first (within 
a factor of 2-3 or so), which puts a limit on the amount of 
ejecta that the fir st generation must p roduc e. Early versions 
of bo th models (|D'Antona fc Caloil l200j ; iDecressin et al.l 
120071 ) suggested non-standard IMFs, heavily weighted to- 
wards the polluter in question, but these IMFs are difficult to 
justify. Current versions suggest that in fact the first gener- 
ation needed to be significantly larger than what is observed 
today, and that the cluster needs to have lost approximately 
90% of its first generation, whil e retaining all of the gas and 
the second generation of stars (Ventura fc D'Antona|[2009l ; 
IDecressin. Baumgardt. fc Kroupal 20081 ) . Both models also 
require a certain amount of dilution of the polluted ejecta 
with primordial gas, in order to match the observed abun- 
dances of light elements, including lithium. 



The work of iPietrinferni et all (|2009l 'l looked at the ef- 
fect of both enhanced helium abundance, and enhanced 
C+N+O abundance, on isochrones appropriate for globu- 
lar clusters. They confirm that multiple main sequences are 
caused by a change in helium abundance at constant metal 
abundance (both iron and light elements). Multiple subgiant 
branches, on the other hand, are caused by a difference in 
total C+N+O ab undance at consta nt helium, constant iron, 
and constant age. lHan et all l|2009T ) show that the split sub- 
giant and giant branches of NGC 1851 are best explained in 
both U — I and V — I by a small change in both helium and 
metal abundance. The alternative explanation is a difference 
in age of the two populations of ~1 Gyr and constant abun- 
dances. It is, however, difficult to explain a delay in star 
formation of such a long time, and difficult to understand 
where the gas for the second generation came from. 



Only a few groups have tried to put together the 
entire scenario . The most successful is the work of 
iD'Ercole et al.l l|2008l ). which incorporated hydrodynamic 
simulations to study the flow of the AGB and super- 
novae ejecta in the vicinity of the proto-cluster. They 
find that the ejecta collects in a cooling flow and re- 
turns to the core of the cluster. They also look at a sce- 
nario in which pristine gas, which was pushed out of the 
cluster vicinity by the supernovae and massive stars, re- 
turns to the cluster after a few million years and mixes 
with t h e ejecta. Similar t o Prantzos, Charbonnel, & Iliadis 



2007); iD'Antona fc Ventural (|2007T l; IVentura fc D'Antona 
2009), they find that the combination of pristine gas 



and ejecta is necessary. They also find that much of the 
first generation in the cluster ma y be lost if th e clus - 
ter i s tidally limited. Finally, both ID'Ercole et al.1 (|2008r ) 
and IDecressin. Baumgardt. fc Kroupal ~ 20081 ) use N-bodv 
dynamical models to study the subsequent evolution of 
this two-generation cluster and the mixing of the two 
populations. They co ncur with the general results of 
iDowning fc Smsl (feoOTh that, at the current time, a two- 
generation cluster will not be dramatically different from a 
single-generation cluster in terms of the dynamics and dis- 
tribution of stars. 



2 A TOY MODEL CALCULATION 

In this paper, we investigate a modified version of the sce- 
narios presented in the introduction. First, we will consider 
the case where both fast-rotating massive stars and asymp- 
totic giant branch stars contribute to the pollution of ma- 
terial which forms the subsequent generation(s). Second, we 
include the possible effects of stellar collisions. 

Recall that the multiple populations phenomenon is 
only seen in globular clusters and not in the field or in open 
clusters. Dense globular clusters are one of the few places in 
the u niverse where stell ar collisions are expected to be preva- 
lent (jHills fc Davlll976l ). There are three relevant outcomes 
of stellar collisions: runaway collisions, modification of the 
rotational properties of stars, and creation of intermediate 
mass stars. 

Dynamical calculations of dense stellar systems have 
shown that there can be dozens of direct stellar collisions 
between massive stars within the first few mi llion years of 
the cluster's life ()Portegies Zwart et al ] |2004T ). These run- 
away collisions involved stars with a total mass of up to 1000 
M© and occur in any cluster which undergoes a core collapse 
withi n the first 3 Myr or so (jFreitag. Rasio. fc Baumgardtl 
2006) The runaway collision dramatically modifies the up- 
per end of the initial mass function. The most massive stars 
in the cluster will not follow their normal isolated evolution, 
but will merge to create the collision product. Therefore, 
there will be fewer supernovae and fewer neutron stars and 
black holes created. Secondly, the evolution of these very 
massive stars must include substantial mass loss, but the 
mass loss history of the collision product is very different 
from that of the sum of the stars which go into the colli- 
sion. As shown bv lGlebbeek et al.l |2009), the total amount 
of mass ejected by the collision product can be up to a fac- 
tor of ~1.5 larger than that ejected by the collision parents 
individually. 

At the same time that the runaway collision is hap- 
pening, other collisions will also be going on in the clus- 
ter. Some of them will produce massive collision prod- 
ucts, and others will be between very small and very large 
stars, which will produce rapidly-rotating stars. Unless two 
stars collide exactly head-on, most of angular momentum 
of the original trajectory is deposited into the collision 
product, producing an object wh ich is spinning rapidly 
l|Lombardi. Rasio. fc Shapirdll996r i. Many of them could be 
rotating with a substantial fraction of their break-up veloc- 
ity. Therefore, the population of fast-rotating massive stars 
should be larger in a cluster with a high collision rate. 

And finally, collisions between low and intermediate 
mass stars could increase the number of stars in the 3-10 
Mq range. Some of the original stars in this range will have 
been removed due to collisions themselves, and so determin- 
ing the total number of intermediate mass stars is not a 
trivial calculation. However, these stars will not segregate 
quickly enough to become part of the sub-cluster which cre- 
ates the runaway collision. This population does not start 
contributing to the gas of the subsequent generation until 
the stars reach the AGB phase, around 30-100 Myr after 
the formation of the cluster. 

To try to get a handle on the abundances of the gas 
which could possibly form the second and third generations 
of stars, we calculate a very simple chemical evolution model 
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based on the scenario described above. We also follow the 
mass function of stars as they form, evolve to remnants, 
and are modified by collisions. We begin with a gasless star 
cluster at t=0. We assume some typical quantities for a glob- 
ular cluster: the total mass of stars in the system is 1 x 10 5 
Mq, its half-mass radius rn=3 pc, and the stellar veloc- 
it y dispersion is 10 km s -1 . Accordin g to the simulations 
of iFreitae. Rasio. fc Baumgardtl l|200rJ l. a runaway collision 
will occur if the cluster has an initial Wo (the dimensionless 
central potential of a King model) of 8 or higher, corre- 
sponding to a concentration parame ter c = log(rt/r c ) = 2 
or higher (|Binnev fc Tremaine!l2008l . figure 4.10). For King 
models, rh /rt is approxi mately constant, with a value of 0.12 
at c = 2 l|Spitzerlll987l . p. 16), and therefore this globular 
cluster must have a core radius r c = 0.25 pc, tidal radius 
r t = 25 pc, and a central density of 2.5 xlO 5 Mg/pc 3 in 
order for this scenario to be valid. These quantities are not 
unreasonable for dense globular clusters seen today. We as- 
sume the stars are formed with a Salpeter mass function 
between 0.1 and 120 M Q and have Z=0.001, Y=0.24 and 
an alpha-enhanced abundance mixture with [a/Fe]=+0.4. 

Early in the cluster's life, we assume that one runaway 
collision occurs. We remove each of the stars which partic- 
ipate in this collision from the initial mass function. The 
yields from this collision are based on the low-meta llicity 
sequence in the calculations of Gl ebbeek et aTl £2009). The 
original calculation used a simplified nuclear network that 
could only follow a small number of species. We have re- 
calculated the yields of the iGlebbeek et al.l (|2009 | '1 calcula- 
tion u sing the nucleosynthesis r outines of Stancliffe et "all 
(l2005h with reaction rates f r om Angulo et al.l ( 19991 ') and 
iHerwig. Austin, fc Lattanzid (|2006h . The reaction network 
includes proton, neutron and a captures on stable isotopes 
up to 34 S and some iron group elements. 

Material leaves the runaway collision products in two 
ways. First, during each collision, the product loses a few 
percent of its total mass because of the energy of impact of 
the collision itself. The amount of each element ejected in 
this way, summed over all collision in the runaway sequence, 
are listed in table[T]as 'ejecta'. Between collisions, the colli- 
sion product is evolving as a relatively normal massive star, 
and exhibits strong mass loss. In our calculations, we assume 
that this loss is from a radiative wind, although LBV-like 
mass loss may also be important. The total amount of mate- 
rial lost by the runaway collision product in winds is listed in 
table [T] There are two wind channels, listed separately: the 
wind of the collision product in between collisions 'Wind'), 
and the wind of the collision product during its remaining 
lifetime after the last collision ('Rem. wind'). These winds 
should leave the parent star with velocities which are typical 
of O stars and luminous blue variables, which range from a 
few hundred to a few thousand km s _1 . In a dense stellar 
system with many other massive stars present, it is plau- 
sible that these winds will collide with the winds of other 
stars, shock, and slow down sufficiently to remain trapped 
in the centre of the cluster and contribute to the second stel- 
lar generation. We assume that is the case initially, but also 
perform the calculation under the assumption that all the 
wind material is lost from the cluster. 

As in Glebbee k et al.l (|2009T l we were only able to fol- 
low the evolution of the collision product about half way 
through core helium burning due to numerical difficulties 



Table 1. Abundances, in mass fraction, from the runaway colli- 
sion. The total amount of mass in each component is given. 



Element 


Eject a 


VV lllQ 


Rem. wind 


Rem. remaining 




1 A x. 1 A/T 
14D.1 IVIq 


ZZ0.4 JVJ.0 


o4.1 IViQ 




X H 


5.97c-01 


3.09e-01 


1.21c-01 


5.70e-02 


4 He 


4.02e-01 


6.90e-01 


8.78e-01 


9.50e-01 


12 C 


1.89e-05 


1.02e-05 


8.78e-06 


7.48e-06 


13 C 


2.16e-06 


2.24e-06 


2.31e-06 


2.12e-06 


14 N 


6.66e-04 


7.10c-04 


7.09e-04 


7.21e-04 


16q 


5.03e-05 


1.79e-05 


1.21e-05 


7.40e-06 


19p 


2.09e-09 


4.44e-10 


1.56e-10 


2.12c-ll 


20 Ne 


7.03e-05 


4.31e-05 


2.85e-05 


2.41e-05 


23 Na 


2.20e-05 


1.29e-05 


3.94e-06 


2.16e-06 


24 Mg 


3.60e-05 


7.50e-05 


9.75e-05 


9.97e-05 


25 Mg 


3.82e-07 


7.65e-08 


3.14e-08 


3.31e-09 


26 Mg 


2.31e-06 


7.35e-07 


2.37e-07 


1.36e-07 


26 A1 


1.95e-06 


1.20e-06 


7.16e-07 


5.73e-07 


27 A1 


4.01e-06 


5.62e-06 


6.11e-06 


6.24e-06 



during this evolution phase. We calculate t he final remnant 
mass o f the collision product as described in lGlebbeek et all 
(2009) and estimate the yields for the remaining core helium 
burning lifetime by assuming that the remainder of the en- 
velope is ejected without undergoing additional processing. 
The results are listed in table [1] under 'Rem. remaining'. 

This particular runaway is only one possible combina- 
tion of parent stars, impact ve l ocitie s, etc. It was chosen 
in the study of IGlebbeek et aTl (|2009h to be representative 
of the collisions seen in N-body simulations. The details of 
exactly how much mass, and the exact composition of the 
ejecta, will change slightly if the details of the runaway col- 
lision changes. This must be kept in mind when interpret- 
ing the results of the chemical evolution calculations pre- 
sented below. A com parison of the t h ree ru naway collision 
sequences studied in iGlebbeek etUI (|2009h show that the 
yields are consistent to within 10-20% depending on the ele- 
ment. The total amount of mass which can be released varies 
by up to a factor of two (up to approximately 1000 Mq). As 
we will discuss below, the results do not change dramatically 
if we adopt the upper limit of these values. 

While the runaway collision is going on, other stellar 
collisions will s pin up some of the massive star s in the clus- 
ter. Following iGiirkan. Freitag. fc Rasid |2004l ). we assume 
that the massive stars in the cluster segregate towards the 
centre, and form a decoupled dynamical cluster which un- 
dergoes core collapse. This cluster, until very close to the 
time of core collapse, has the same core radius and central 
density as the original cluster. However, the mass function 
becomes much more weighted towards massive stars, so that 
the average stellar mass in this sub-cluster is more like 20 
Mq (compared to ~ 0.35 Mq for a normal cluster). Us- 
ing these parameters a nd the equation for the average time 
between collisions from lLeonar d (1989), we find that approx- 
imately 1100 stellar collisions should occur during this first 
5 Myr in this initial sub-cluster of stars. There are approx- 
imately 200 stars in this cluster with masses above 20 Mq, 
forming approximately 35% of the stars in the sub-cluster. 
Therefore, it is likely that every high-mass star has under- 
gone at least one stellar collision. 

Following iDecressin et all l|2007h . we define "rapidly ro- 
tating" to mean having a rotation rate that is at least 80% 
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of the critical rotation rate. We expect that some fraction 
of massive stars are primordial fast rotators (i.e. they were 
b orn that way and not s pun up by collisions). The models 
of iDecressin et all l|2007f ) assume that all massive stars are 
rotating rapidly enough to contribute polluted material to 
the secon d generation. Obse rvations of stars in young clus- 
ters (e.g. iDufton et al-llioO^ 1 ) suggest that cluster stars are 
more rapidly rotating than those in the field, and that the 
fraction of rapid rotators is more like 20%. Initially, we make 
the extreme assumption that every collision will turn a slow 
rotator into a rapid rotator. This assumption is unlikely to 
be completely correct, as the amount of angular momentum 
which can be added to a high-mass star depends on the mass 
of the impactor, its velocity and its position of impact on the 
high mass star. However, the highly collisional environment 
of the mass-segregated sub-cluster means that high mass 
stars will likely undergo more than one collision. The net 
effect of these collisions will be to spin up the population. 

We use the yields of lDecressin et all l|2007h to determine 
the contribution of this population to the ejected gas. We 
use the yields of their case C reaction rates for 40, 60 and 
120 M© fast-rotating stars. These yields come from a set 
of reactions in which the reactions involving 20 Ne through 
27 A1 are set to their experimental upper or lower limits in 
such a way to produce the most favourable set of yields for 
this work. However, we do not use the yields from the mod- 
els in which the 24 Mg(p,7) reaction rate is increased by a 
factor of 1000 at 5 x 10 7 K (case D). This fudge was re- 
quired to produce the lower range of the observed [Mg/Al] 
values in NGC 6752. However, for t his toy model, w e prefe r 
to use the 'standard' set given by IDecressin et aU (120071) . 
The models presented in that paper assumed a metallicity 
of Z=0.0005 and are initially alpha-enhanced. The yields for 
the runaway collisions and the AGB stars are for Z=0.001, so 
we should be using different yields to be consistent. Unfor- 
tunately higher metallicity yields for fast-rotating massive 
stars are not available. We do not expect our results to be 
dramatically different if we had used other reasonable yields 
or models. 

Like the runaway collisions, these fast-rotating massive 
stars also lose material in two ways. One, a slow wind is 
present, primarily from the equator of the stars, as the sur- 
face material reaches an angular velocity higher than the 
local stellar escape velocity. This material will remain in the 
cluster, and those yields are assumed to be those up to the 
time t he star has reached the end of central H burning, fol- 
lowing [Decressin~e^al. (2007). These stars also have a fast 
wind, with velocities typical for O and B stars (a few hun- 
dred km s _1 ), and is dominant between the end of central 
H burning and the time when the helium-burning products 
appear at the surface of the fast-rotating stars. We make 
the same assumptions as for the winds of the runaway col- 
lision: initially we assume these winds shock and remain in 
the cluster, and then we will present a calculation in which 
they are assumed to be removed entirely from the cluster. 

The time of the first supernova is set by the most mas- 
sive star which is not involved in a runaway collision, or the 
lifetime of the runaway itself, whichever comes first. At this 
time, which is approximately 5 Myr, we assume that all the 
gas in the cluster forms the second generation of stars, and 
any subsequent gas is removed from the cluster until the 
AGB stars begin to contribute. For this reason, we neglect 



the contribution of rapidly rotating 20 Mq stars, as their 
main sequence lifetimes are closer to 10 Myr. We calculate 
the total mass in ejecta from both the runaway collision and 
the fast rotating massive stars, under the assumption that 
both fast and slow winds are retained in the cluster. We find 
that we have 3459 Mq of material, which we assume forms 
the second generation of stars with 100% efficiency. We pop- 
ulate a Salpeter IMF from 0.1 to 120 Mq, as was done with 
the initial cluster. 

The abundances of this second generation are quite ex- 
treme. In Figure [TJ we show the helium abundance of each 
generation, labelled as "primordial", "runaway & FRMS" 
for this second generation, and "all AGB" and "high mass 
AGB" for different possibilities of the third generation, to be 
discussed below. This figure also shows the calculated num- 
ber of low mass stars, with masses less than 0.8 Mq, which 
are expected to be observable members of the cluster today. 
Note that we plot the logarithm of the number of stars. The 
runaway collision, which contributes only 539 Mq of ma- 
terial to the cluster, is strongly enhanced in helium, with 
an overall Y of 0.68. The fast-rotating massive stars are 
also strongly helium enhanced, with Y=0.43 for the 2920 
Mq ejected. This results in an overall helium abundance 
for this first generation of Y=0.47. This is much higher than 
is inferred in even the most extreme globular cluster second 
generation. 

The abundances of other elements, however, match the 
observations reasonably well. Figures [2] and [3] compare the 
calculated abundances of each generation or pollution mech- 
anism, as labelled, to the observations. The shaded region on 
each diagram encompasses most of t he observed ab u ndance 
tren ds, as taken from the f igures in ICarretta et al.l (2009a) 
and ICarretta etakl (|2009bh . The runaway collisions alone 
seem to be the only way to reach significant oxygen deple- 
tion ([O/Fe] ~ -0.9) and sodium enhancement ([Na/Fe] ~ 
1.0). The fast- rotating massive stars do deplete oxygen and 
enhance sodium as well, although to a lesser extent. Because 
the runaway collision produces so little mass, the overall 
abundances of the second generation are oxygen-poor and 
sodium-rich at levels which are consistent with all but the 
most extreme observations in clusters (e.g. in NGC 2808). 

Magnesium and aluminum abundances from this second 
generation are also consistent with most of the observations. 
Again, none of our models are able to reproduce the most 
extreme populations (the high aluminum, low magnesium 
portion of the diagram). It should be noted, however, that 
there are very few stars in that region and they are all in 
NGC 2808; the bulk of the stars in other clusters lie in the 
vertical region of this plane, and we argue that magnesium 
depletion needs to be confirmed in more clusters before this 
diagnostic is used to rule out particular models. We also note 
that this generation does not produce very much aluminum, 
despite the interestingly high helium yields and the good 
agreement with the bulk of the Na-0 observations. 

Next, we consider the likelihood that collisions in this 
early stage will produce additional intermediate-mass AGB 
stars later in the cluster's life. We used the properties of the 
entire cluster core rather than that of the sub-cluster of mas- 
sive stars because the low to intermediate mass stars that we 
were considering her e have no t had t ime to mass segregate. 
Using equation 13 of iLeonardl (|l989h . we found that only 4 
collisions should have occurred in the first 5 Myr. Even if 
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we assume that the cluster conditions remain the same until 
the stars of this mass reach the AGB phase (about 30 Myr), 
only a few dozen collisions will occur, and only approxi- 
mately 5% of those will produce stars with masses between 
3 and 8 Mq. Therefore, we will neglect the contribution of 
collisionally-created AGB stars. 

Now we calculate the contribution of the AGB stars 
from both the first and second stellar generations. We use 
yields for Z =0.001 for both the first and seco nd generation 
populations (|Ventura fc D'Antonall200sl . l2009l ). These yields 
will be incorrect for the AGB stars formed from the ejecta 
of the runaway collision and fast rotating massive stars, as 
the helium and light element content of these stars is ini- 
tially increased. The structure of helium-rich stars is dif- 
ferent than that of normal stars, and the hydrogen-burning 
cycles are dependent on the abundances of the catalyst ele- 
ments. Therefore, we should not simply scale the yields from 
the normal AGB stars to estimate the yields of high-helium 
stars. The helium abundances will be higher than given by 
the normal star yields because even unprocessed material 
will be enriched in helium. However, the number of second- 
generation AGB stars is low (~ 65), and so we will take 
the conservative assumption that their yields are the same 
as the first generation population. In addition, helium-rich 
AGB stars have lifetimes which are ~ 50-70% shorter than 
stars of the same mass but normal helium. In our simple 
model, this does not affect our results because we simply 
sum up all the contributions from AGB stars with lifetimes 
shorter than our cutoff. However, a more detailed model will 
need to include this lifetime effect. 

First, we assume that all AGB stars between 3 and 6 
Mq contribute to the material which forms the third genera- 
tion. The AGB stars produce almost twice as much material 
as the massive stars (6084 Mq). The helium abundance of 
this population is Y=0.29, which is higher than the stan- 
dard value but still not as high as the Y=0.4 inferred for 
clusters such as NGC 2808. The AGB ejecta is slightly en- 
hanced in oxygen compared to the initial value ([O/Fe] = 
0.55, up from 0.4), and is significantly enhanced in sodium 
([Na/Fe] = 1.04, up from -0.2). The sodium enhancement is 
much larger than that seen in most globular clusters, and the 
AGB stars alone do not produce stars with low oxygen/high 
sodium values that are seen in globular clusters. Dilution 
with primordial material reduce both the sodium and oxy- 
gen abundances, but an oxygen depletion of approximately 
1 dex is impossible to accomplish with these AGB yields. 
Similarly, this population produces some aluminum without 
much change in magnesium. This population is labelled 'all 
AGB' in figures [T] -El 

However, 3 M AGB stars have lifetimes of over 300 
Myr. This time is long enough that SNIa may have started 
to pollute the cluster and disrupt the gas. Also, the most 
massive AGB stars started losing their mass after only ~ 
50 Myr, and it is not clear that this material would have 
remained in the cluster, waiting for the ejecta of the lower 
mass stars. It is more likely that the longest possible time 
for the AGB ejecta to collect is more like 100 Myr, which is 
the lifetime of a 5 MqAGB star. If we restrict ourselves to 
only the most massive AGB stars (5-6 M ), then the sodium 
and oxygen yields are more consistent with the observations, 
and in fact are very similar to those from the runaway + 
FRMS population. Under this assumption, AGB stars only 



contribute ~ 2100 Mq to the new generation, an amount 
of mass which is comparable to that of the first generation. 
This population also produces more aluminum and shows a 
very slight magnesium dilution. This population is labelled 
'high mass AGB' in figures Q] - 

If we combine the second and third generation, we have 
9.5 x 10 3 Mq of material, or less than 10% of the initial 
cluster. If we compare the number of stars that will still be 
in the globular cluster at current time (less massive than 
0.8 Mq), the two new generations have created 25 000 stars 
and there are approximately 2.6 x 10 stars from the first 
generation. We are still required to lose 90% of the initial 
low-mass stars in order to have our younger generations form 
half the cluster at the current day. These numbers assume 
that all AGB stars contribute to the third generation, which 
is the most generous assumption one can make about total 
mass, but is almost certainly an overestimate as discussed 
above. 

The other way to mitigate this mass problem is to allow 
primordial material to mix into the gas which will form ei- 
ther the second or third generation (or possibly both) . While 
this will certainly help boost the mass of that generation, 
it will also change the abundances. In figures [2] and [3] we 
have drawn lines of dilution for the second and the two pos- 
sible third generations. The amount of dilution ranges from 
almost none near the points labelled by the polluters, to a 
huge amount of mass as the line nears the primordial abun- 
dance. We feel that it is more likely that the second genera- 
tion would be polluted than the third, and so we calculated 
the total amount of mass needed to bring the helium abun- 
dance down to some values of interest. For example, only 
1500 M Q is needed to bring the helium abundance to Y=0.4, 
but 77 000 M brings Y to 0.25. This is comparable to the 
mass of the initial generation. 

The previous results assume that both the "slow" and 
"fast" winds from runaway collision products and fast- 
rotating massive stars contribute to the second generation 
in the globular cluster. Fast winds, with velocities between 
a few hundred and a thousand km/s, may not be retained in 
the potential well of the initial cluster. We have calculated 
the total mass of gas, and its abundances, using only the slow 
winds from the two contributors to the second generation. 
As expected, the total amount of mass available is reduced 
(« 1900 Mq instead of « 3500 Mq), and the abundances 
are different. In particular, the helium content is much lower 
(Y=0.321) because the fast wind is thought to be occur- 
ring when the stars are predominantly helium stars. The 
oxygen and sodium abundances are both approximately 0.1 
dex lower, while the aluminum and magnesium abundances 
are lower by more like 0.15 dex. These abundances are still 
consistent with the observed abundances in clusters. The 
oxygen and sodium abundances are closer to the extreme 
values observed, while the aluminum/magnesium values are 
closer to the primordial values. 

For comparison, we have also calculated the expected 
contributions from our three pollution processes if we ne- 
glect the effects of stellar collisions entirely. In this case, 
runaway collisions do not contribute at all. We assume that 
only 20% of massive stars are rotating rapidly enough to 
produce any polluting material (but we do include both 
slow and fast winds from fast-rotating massive stars). We 
assume that the other 80% of massive stars are not rotating 
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Figure 1. Histogram of helium abundance of the three stellar 
generations, and the number of low-mass (i.e. observable) stars 
in each generation. The polluter(s) of each generation is marked: 
primordial abundance, all asymptotic giant branch stars, high 
mass asymptotic giant branch stars only, and runaway collision 
product & fast-rotating massive stars. Note that the y axis is 
logarithmic. 



Figure 3. Aluminum vs. Magnesium abundances for the three 
stellar generations. The sh aded region shows the extent of the 
observations, as taken from [Caxretta e t al. (2009a b). The abun- 
dances of the runaway collision products, the population of fast 
rotating massive stars, and the population of high-mass AGB 
stars only, are also marked. The curved lines connecting the gen- 
erations to the primordial abundances show the effects of dilution. 
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Figure 2. Sodium vs. Oxygen abundances for the three stellar 
generations. The sha ded region shows the e xtent of the observa- 
tions, as taken from [Carrctta et al. (200930)- The abundances 
of the runaway collision products, the population of fast rotating 
massive stars, and the population of high-mass AGB stars only, 
are also marked. The curved lines connecting the generations to 
the primordial abundances show the effects of dilution. 



at all, and we neglect any contribution of their winds to the 
second generation. This is reasonable because their surface 
abundances have not been modified by nuclear burning. The 
effects of AGB stars are reduced as well, since we have much 
less mass with which to make second-generation AGB stars. 
Therefore, the total amount of mass in the subsequent gen- 
erations is only 6500 M©, a mere 6.5% of the initial mass 
of the cluster. Of this material, the bulk is from the AGB 
stars of all masses, and so is very high in oxygen, sodium, 
and aluminum, and only mildly enhanced in helium. 



3 DISCUSSION AND CONCLUSIONS 

In this paper, we take another theoretical look at the prob- 
lem of multiple populations in globular clusters. In particu- 
lar, we tried to address the observational constraint that this 
is specifically a dense old cluster phenomenon. We concen- 
trated on the effects of stellar collisions, which is one of the 
only physical mechanisms which affects clusters much more 
strongly than anywhere else in the universe. We also com- 
bined two previous models for multiple populations - the 
fast-rotating massive stars and the asymptotic giant branch 
stars - to determine if having more polluters would help with 
either abundances or with the total mass in the subsequent 
generations. 

We included the effects of collisions in three ways: the 
impact of a runaway collision, an increase in the number of 
fast-rotating massive stars, and an increase in the number 
of intermediate mass AGB stars. We also allowed for two 
subsequent generations. One is formed within a few million 
years of the first cluster formation, from the ejecta of run- 
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away collisions and fast-rotating massive stars. The other 
is formed approximately 100 million years later, from the 
ejecta of the AGB stars. 

Our results are not substantially different from those of 
other groups who are modeling the early evolution of multi- 
ple population clusters. We still find that the mass produced 
by the polluters alone is insufficient. Loss of a large fraction 
of the initial cluster, dilution by primordial gas, or both, 
is required. However, stellar collisions do produce material 
with abundances which are interestingly extreme in all el- 
ements, and may well be a substantial contribution to the 
extreme populations in the most massive clusters. 

It is clear, from this work and others, that we do not yet 
have a complete understanding of the multiple populations 
issue. Here, we try to highlight a number of the largest re- 
maining problems, starting with those specific to our model 
and then discussing the puzzle more generally. 

So far, there have been no simulations of star formation 
in the environment (s) envisioned here. Star formation typ- 
ically occurs in cores in molecular clouds, an environment 
which is cold, dark, dense, and shielded from the outside 
universe. In the centre of a dense cluster, a cloud of ejecta 
will be subject to the radiation field of the first generation 
stars. No simulations of star formation in a region with this 
kind of external radiation bath have been done, and so we 
do not really understand the properties of the second gen- 
eration. We have assumed a Salpeter mass function for our 
subsequent generations, for example, with a maximum mass 
of 120 Mq . It may be reasonable to postulate that high mass 
stars coul d never form in the sec ond generation (as was sug- 
gested by iD'Ercole et al l (|2008l )). which would mean that 
the second generations have more low-mass stars then in 
our toy model. For a Salpeter mass function between 0.1 
and 120 Mq, 44% of the mass is found in stars more mas- 
sive than 0.8 Mq. Under the assumption that only stars of 
mass 0.8 Mq or less form in the subsequent generations, 
we can increase the predicted mass of these generations by 
approximately a factor of two. 

If the cluster needs to be substantially more massive 
initially than assumed in this paper, then we should consider 
how the different polluter channels scale with cluster mass. 
For fast-rotating stars and AGB stars it is reasonable to 
assume that the yields scale with the total cluster mass. 
For the collision runaway this is not so clear. The mass of 
the collisional runaway as a function of the cluster mass 
has not been considered in detail in the literature. If the 
mass of the collision runaway scales with cluster mass, then 
it behaves as the other two polluter models. However, if 
the mass of the collision runaway depends more steeply on 
the mass of the cluster, then it will become relatively more 
important with increasing initial cluster mass. In that case 
the number of stars that need to be lost from the cluster can 
be smaller than when either of the other two scenarios alone 
are considered. If the total mass of the runaway is only a 
weak function of cluster mass the other polluter scenarios 
will become relatively more important for increasing cluster 
mass. 

In these models, we have neglected the effects of bi- 
nary stars. We know that globular clusters do have bi- 
naries, although work suggests that the fraction may be 
lower than in the field and open clusters (|Davis et al ] |2008l ; 
iRubenstein fc Bailvnl [l997h . Binary stars can modify our 



toy model in a number of ways. Collisions betwee n binary 
stars are more likely than between single stars l|Leonardl 
Il989h because of the larger cross section of the binary, and 
many of those interactions can result in more than two stars 
merging. Binaries can also incre ase the likelihood of having 
more than one runaway collision (Gur kan. Fregeau. fc Rasiq 
2006). There are indications that massive stars h ave a higher 
binary fraction than low mass stars in the field l|Ladal 2006) 
and so it may not be unreasonable to expect that interacting 
binary stars may have a significant impact o n the mass lost 
from massive stars. Follow ing suggestions bv lde Mink et all 
l|2009h , IVanb evcrc u (2009) goes as far to suggest that inter- 
acting massive binaries are responsible for all the pollution 
in globular clusters, not fast-rotating massive stars or AGB 
stars. 

All models of multiple populations need to be refined 
substantially in order to explain the cluster-to-cluster varia- 
tions. While it is clear that something is going on in almost 
every well-studied cluster, it is not clear that we understand 
how that effect depends on the cluster properties. In our sce- 
nario, we would argue that the cluster-to-cluster variations 
were caused by different collision rates early in the clusters' 
lives, perhaps driven by slightly different formation condi- 
tions and initial densities. In fact, it may be that we can use 
the extreme abundances produced in runaway collisions to 
determine which clusters hosted a runaway all those years 
ago. Only some clusters show evidence for extremely high 
Y values (Y ~ 0.4) and very extreme sodium-oxygen anti- 
correlations, such as ui Centauri and NGC 2808. Perhaps 
only those clusters were formed with sufficiently high initial 
density to induce a runaway. A more detailed study of the 
abundance patterns in the individual pollution mechanisms 
may help disentangle these processes. 

It has become clear that the epochs of globular clus- 
ter formation and very early evolution are crucial pieces of 
the puzzle when disentangling the present-day properties of 
these ancient objects. It is also clear that our understanding 
of the dominant processes and effects during these epochs is 
not as strong as we would like. Learning about the early life- 
times of dense clusters is not easy - the stellar archaeology 
required is quite intricate. We are looking at small changes 
in current surface abundances and brightnesses of very old, 
low mass stars, and inferring a significant amount of action 
involving more massive stars over 10 billion years ago. How- 
ever, we have learned a great deal about this phenomenon 
since it first was identified only a few years ago, and progress 
will continue to be made. 
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